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ABSTRACT 

We establish a classification scheme for stochastic acceleration models involv- 
ing low-frequency plasma turbulence in a strongly magnetized plasma. This clas- 
sification takes into account both the properties of the accelerating electromag- 
netic field, and the nature of the transport of charged particles in the acceleration 
region. We group the acceleration processes as either resonant, non- resonant or 
resonant-broadened, depending on whether the particle motion is free-streaming 
along the magnetic field, diffusive or a combination of the two. Stochastic accel- 
eration by moving magnetic mirrors and adiabatic compressions are addressed 
as illustrative examples. We obtain expressions for the momentum-dependent 
diffusion coefficient D(p), both for general forms of the accelerating force and 
for the situation when the electromagnetic force is wave-like, with a specified 
dispersion relation u = u(k). Finally, for models considered, we calculate the 
energy-dependent acceleration time, a quantity that can be directly compared 
with observations of the time profile of the radiation field produced by the accel- 
erated particles, such as during solar flares. 

Subject headings: Sun: activity - Sun: flares - Sun: particle emission - Sun: 
X-rays, gamma rays 

1. Introduction 

The acceleration of particles to non-thermal energies is a fundamental phenomenon 
throughout the Universe. In particular, the impulsive phase of a solar flare is characterized by 
the acceleration of charged particles, including electrons, protons and heavier ion species (see, 
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At the fundamental level, of course, charged particles are accelerated by an electric field; it 
is the degree of spatio-temporal coherence of this electric field that distinguishes acceleration 
mechanisms. 

The action of an accelerating field can be divided into two main categories, depend- 
ing on whether it relies on systematic or stochastic change in the energy of the particles. 
Systematic change of energy is involved in a a broad range of particle acceleration mecha- 
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On the other hand, when the energy change has a fluctuating character, the mechanism 
of particle a cceleration i s said to be stochastic. The classic example is second-order Fermi 



acceleration (jFermilll949l ). which is based on a more random motion of the scattering centers 
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Extended w ork on particle acceleration in solar flares has shown that stochastic accelera- 
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not only avoid some of the fundamental large-scale electrodynamic issues characteris- 



tic of systemic acceleration models (lEmslie k Henoux 



19951 : Miller et al.lll997f ). but also can 



produce an acceleration efficiency (lEmslie et al.l 120081 ) that is broadly consistent with that 
deduced from hard X-ray observations. Given, then, the likelihood that stochastic accelera- 
tion plays an important role in the acceleration of non-thermal particles in solar flares, we 
seek in this work to synthesize various stochastic acceleration models in a unified manner. 
In this way, we highlight their similarities and differences, and we show how their essential 
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paramet ers could be constrained through observat ion of, e.g., turbulent and di rected plasma 



motions (lAntonucci et al.lll982t iFludra et al.lll989f ) or of magnetic fluctuations (IKontar et al. 
2011bl ; iBian et al.ll201ll ) in solar flares. 



We classify stochastic acceleration models by the nature of the spatial transport of 
particles in the acceleration region. We do this through the ICorrsinl ( 119591 ) approxima- 
tion, a procedure for relating Eulerian and Lagrangian corre lation functions which enjoys 



widespread use in the theory of turbulent transport (s ee, e.g., IWang et al.lll995t IVlad et al. 



2004 ; iNeuer fc Spatschekl 120061 ; iTautz fc Shalchil |2010| ) . This approach allows us to divide 
stochastic acceleration models into three primary classes: resonant, resonant-broadening and 
non-resonant, depending on the nature of the transport process affecting the accelerated par- 
ticles. The necessary analysis is presented both through an ordinary space representation 
and in the Fourier domain, a representation particularly appropriate when the acceleration 
mechanism involves waves. For each case we derive the form of the diffusion coefficient D(p) 
(when it exists) as a function of the particle momentum p. We then use this to obtain 
the characteristic acceleration time as a function of particle energy E in the non-relativistic 
regime. 

We also present examples of scenarios that depend on the specific nature of the force 
accelerating the particles, with the turbulent fluctuations being associated either with weak 
wave turbulence or strong turbulence. It is emphasized that stochastic acceleration can also 
occur when the accelerating field is coherent, provided the particle transport in this field is 
itself stochastic, for instance diffusive. 

In Section [21 we revie w the fundament al concepts of stochastic acceleration mechanisms 
through a discussion of the lObukhovl (119591 ) model, one of the first mathematical formulations 
of turbulent acceleration. In Section [3, we derive the essential features of resonant, non- 
resonant, and resonant-broadened models, respectively. Examples of different forms of the 
accelerating force are considered in Section HI A brief discussion of the energy spectra that 
result is given in Section 5 and an overall summary is presented in Section [6j 



2. Stochastic Acceleration 
2.1. Fundamental Concepts 

2.1.1. Velocity and Spatial Diffusion 

Consider a particle of momentum p subjected to a rapidly fluctuating force F(t), so that 
the (one-dimensional) motion of the particle is described by a set of Langevin equations 
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p = F(t) ; mx — p = mv . (2-1) 

The random force F(t) is taken to be a Gaussian process with zero average < F(t) > = 
and non-zero variance < F 2 >. The force F(t) is delta-correlated in time: 



C(t, f) = < F(t) F(t') > = 2D 5(t - t') , (2.2) 

where D is a diffusion coefficient and the averaging is made over ensemble realizations. 
As a consequence of stationarity, the (auto-) correlation function has the property that 
C(t,t') =C(t-t'). 



The lObukhovl (119591 ) model (12. ip is the undamped version of the Ornstein-Uhlenbeck 
process, described by 



p = F(t) ; mx = p , (2.3) 

r 

where r is a damping time. From the corresponding Fokker-Planck equation for the dis- 
tribution P(p,t) of momentum at time t, one finds that the stationary distribution has a 
Gaussian form 



P(P)= /^^= exp f-^W) > ( 2 - 4 ) 
a/2tt<j9 2 > V 2<p 2 >J 

where < p 2 >= tD is the variance characterizing the momentum distribution. The mean 
square displacement grows linearly with time, < x 2 >= 2D x t, where the spatial diffusion 
coefficient D x is related to the momentum damping rate by D x = r(< p 2 > /m 2 ). Spatial 
diffusion occurs only for thermalized particles when the statistics of the process have become 
stationary, i.e., when < p(t)p(t') > depends only on (t — if) according to < p(t)p(t') > = < 
p 2 > exp(— | t — t' | /t). 

On the other hand, the undamped (r — > oo) Obukhov model results instead in mo- 
mentum diffusion: the mean energy of the accelerated particles, < p 2 >, grows linearly 
with time. The equation describing the evolution of the probability density P(x,p,t) =< 
S[x — x(t)] 5[p — pit)} > is the Fokker-Planck equation 



dP(x,p, t) 



p dP(x,p,t) 
m dx 



D 



d 2 P(x jPl t) 



(2.5) 



dt m dx dp 2 

With the initial condition P(x,p, t = 0) = S(x) S(p), the solution of this equation is a Gaus- 
sian form exhibiting simultaneous unbounded spreadi ng with respect to both m omentum 
and position and with the following scaling behaviors ( [Thomson &: Benfordlll972[ ): 
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2Dt 3 

< p 2 >= 2Dt ; < x 2 >= —- . (2.6) 

6m z 

Averaged over position x, the probability density P(p, t) is a Gaussian with a variance that 
grows linearly with time. Similarly, averaged over momentum p, the probability density 
P(x, t) is a Gaussian with a variance that has a cubic dependence on time^. 

It should be emphasized that the standard diffusion law < p 2 > oc t produced by this 
model is a particular case. In general, stochastic acceleration processes result in a momentum 
diffusion characterized by < p 2 > oc t a , with a / 1. Cases with < a < 1 are termed 
sub-diffusive and cases with a > 1 are termed super- diffusive (jBalescul 120001 ) . Anomalous 
diffusion occurs when the momentum diffusion coefficient is zero, infinite, or if it depends on 
momentum. In general, the momentum diffusion coefficient is momentum-dependent when 
(a) the force itself depends on momentum, or (b) the force depends on position as well as 
on time. 



2.1.2. Spatially-Varying Force Fields and Correlation Times 

We now consider a non-trivial extension of the Obukhov model which incorporates an 
explicit spatial dependence of the force F(x,t): 

p = F(x,t) , mx = p . (2.7) 

The force is still assumed to be homogeneous and stationary, with the statistical property 
< F(x,t) > = 0. Its Eulerian correlation function is 

C(x, t) =< F(0, 0) F(x, t) > . (2.8) 
We next define the Eulerian spatial and temporal scales by 

I poo i poo 

T= cml dtmty ' x= cmL dxC(x < 0) < (2 - 9) 

respectively, where 



1 The dispersions in velocity and in space are related by < p 2 > oc < x 2 > 1 / 3 , a result that can be used to 
dedu ce the K o lmoro gov law E(k) tx fc -5 / 3 of hydr odynami c turb ulence. This was the original motivation for 
the lObukhovl (|l959h model. Later, it was used by iDupred (]1966f ) in his foundation of resonance broadening 
theory applied to wave-particle interactions in plasmas. 
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C(0,0)=<F 2 > . 
A general expression for the Eulerian correlation function is therefore 



(2.10) 



C(x,t)=<F 2 > c(j,~ 
where the non-dimensional function c satisfies the following conditions: 



(2.11: 



£ (0,0) = l; Jj{^)iQ^l; 
A suitable functional satisfying these properties is the Gaussian 



c[0,-)d 



(2.12) 



x t 
' A' r 



(2.13) 



With the spatio-temporal statistics of the force now known, the problem is to calculate 
the momentum diffusion coefficient (w hen the latter exists) and to determine the form of 
the distribution function. iTaylorl (Il922[ ) established that the momentum diffusion coefficient 
D(p) is related to the time integral of the Lagrangian correlation function, viz. 



Dip) 



dtC L (p;t) 



(2.14) 



where the Lagrangian correlation function of the force is defined by 



C L (p;t) =< F(0,0)F(x(t),t)> 



(2.15) 



Taylor's result (EHJ) follows from < p 2 > = jj dt' J* dt" < F(t')F(t") > = 2 f* dt' C L (t') (t - 
t') and letting t — » oo (when the integral converges). The Lagrangian correlation time tl, 
the correlation time associated with the force felt by a particle along its trajectory, is defined 



as 



Cl(0) 



dtC L (p;t) 



(2.16) 



2 Note the appearance of a factor of four in the temporal Gaussian exponent; this is a direct consequence 
of the half-space limits [0, oo) of the temporal integral (versus the full-space limits (—00,00) in the spatial 
integral) in Equations (|2.9[) . 



-7- 



Since < F 2 > = C(0,0) = C L (t = 0), Equations (12TT0D . (12TT4D and (jgJSD show that the 
diffusion coefficient, when it exists, is given by 

D(p) = r L <F 2 > . (2.17) 
The special cases D = (sub-diffusion) and D — > oo (supe r-diffusion), correspondin g to tl 



being zero or infinity, respectively, have been discussed by iBian fc Browning (120081 ) on the 
basis of a generalized Obukhov model. 

In general, both tl and < F 2 > may depend on the particle momentum p, so that 
D = D(p). We shall find that a typical dependence D(p) is a power law 

D{p) = D p a . (2.18) 
Then, since d < p 2 > /dt = 2D(p) = 2Dop a , it follows that 

< p 2 >cx t 2 '^ , (2.19) 

for a < 2. For a > 2, < p 2 > can attain infinite values within a finite period of time as may 
be seen by solving the equation dp 2 /dt = 2D p a with a > 2. 

It is clear from Equation (I2.17P that determination of the diffusion coefficient essentially 
involves determination of the Lagrangian correlation time tl. When the force is independent 
of position (as in the Obukhov model), determination of ti is a trivial task, because in this 
simple case tl = r, the Eulerian correlation time. A much more significant challenge is to 
obtain the Lagrangian correlation scales (as a function of the Eulerian ones) when the force 
depends on position as well as time; indeed, it is this very spatial dependence of the force that 
makes the stochastic acceleration problem inherently non-linear. A general strategy can be 
summarized as follows: we seek the functional dependence of the Lagrangian correlation time 
Ti (and hence D) on the Eulerian quantities r, A and possibly^ < F 2 >. From dimensional 
arguments, if the time scale tl depends on the length scale A this dependence is expected 
to also involve quantities related to the spatial transport of particles, such as the particle 
momentum p (if the particles are free-streaming) or the particle diffusivity coefficient k (if 
the transport of particles involve some form of spatial diffusion). The pertinent transport 
time scales are mX/p and A 2 /k, respectively. 



Dupred (1966) was the first to propose a scheme, based on the Obukhov model, to include the functional 



dependence of tl on < F 2 >. 
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2.1.3. Remarks 

It must be noted that the above analysis is one-dimensional. Nevertheless, it is still 
a valid description of the behavior of a distribution of charged particles moving under the 
influence of random electric fields in a volume permeated by a uniform magnetic field B , if 
the x-coordinate is taken to represent the coordinate parallel to the magnetic field and the 
accelerating forced is expressed as a parallel electric field: F x {t) = qE(x, t). It is well known 
that the presence of a strong background magnetic field permits the decomposition of the 
motion of the charged particles into a guiding-center motion parallel to the x-axis and a 
much faster perpendicular gyration around the guiding magnetic field line. At this stage, we 
identify the position x of the particle with that of the gyrocenter; this approximation will be 
relaxed later (Section 14.21) . The basic mechanisms for perpendicular acceleration of particles 
by tra nsverse electric fields can also be outlined within the context of a one-dimensional 



model (jSturrock!ll966l ) 



In this work we restrict ourselves to the consideration of low-frequency turbulence and 
so do not consider the effect of gyroresonances. Further, we employ a test-particle approach, 
so that effects of the collective behavior of the particles are neglected. The dissipation of 
electromagnetic energy resulting from the stochastic acceleration of particles will be discussed 
in a separate publication. 



2.2. The Corrsin approximation 

Let us write the Lagrangian correlation function (12.151) in the equivalent form 

C^ t) = J dx <F M F (xAS[x -m>, (2.20) 

where x(t) is a solution of the equations of motion (12.71) . One way of obtaining a relation 
between the Lagran gian corr e lation Cl{p~, t) and the Eulerian correlation C(x, t) is to invoke 



a procedure due to ICorrsinl (119591 ). in which x(t) is replaced by its statistical average, so 
that we may replace 8[x — x(t)] in Equation (12.201) by < S[x — x(t)] >. This leads to the 
factorization 

C i( p;^/*<F(0,0)F(M)><^-, W ]> (2.21) 



4 Othcr possible types of force fields are discussed in Section 4. 
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and immediately allows us to write 



C L (p;t)= dxC(x,t)P(x,t) , 



(2.22) 



where P(x,t) =< 8[x — x(t)] > is the probability density for a particle to move from its 
starting point at the origin to position x over time t. Equation (I2.22p shows that the 
Lagrangian correlation function is the spatial integral of the product of two quantities: 

1. C(x,t), the Eulerian correlation function, which depends on the properties of the force 
field F(x, t), and 

2. the probability function P(x,t), describing the spatial transport of particles in the 
acceleration region. 

This separation into properties of the accelerating force and properties of the particle trajec- 
tories forms the basis for our classification of acceleration models. Specifically, we shall cat- 
egorize stochastic acceleration models primarily through the nature of the quantity P(x,t), 
and then sub-categorize them according to the nature of the force F(x,t) [and hence the 
Eulerian correlation C(x,t)] acting on the particles. 

With the Lagrangian correlation function Cx(p;t) defined as in Equation (12.221) . the 
momentum diffusion coefficient D(p) may be written 



At this stage, a preliminary comment can be made. First, let us use Equation (12.111) to write 
(I2.23P in the equivalent form 




(2.23) 



D(p) =<F 2 > 



I 



oo 




P{x,t) 



(2.24) 



We notice that when 



P(x,t) ~ 8(x) , 



(2.25) 



then 




(2.26) 
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This is essentially the Obukhov result t l = r; the spatial dependence of the force, and hence 
the spatial transport of the particles, is irrelevant in evaluating the acceleration efficiency 
because the temporal variation of the force is taken to be infinitely fast, i.e., r — > 0. For such 
a situation, all stochastic acceleration models reduce to the common limit corresponding to 
the Obukhov model. 

However, there are certainly situations where it is unreasonable to assume that the 
Eulerian correlation time r — > 0. An obvious example is when the force field acting on the 
particles is produced by wave motions, for which case the Eulerian correlation time r is of 
the order of the wave period T. This in turn means that the Lagrangian correlation time Tl 
now has a dependence on the Eulerian spatial scale A, because the particles of necessity also 
feel the spatial variation of the force. 

When acceleration results from the action of waves with a well-defined dispersion relation 
uu — oo(k), it can be more convenient to express the momentum diffusion coefficient in terms 
of Fourier components of the Eulerian correlation function of the force field, i.e., in terms of 
the spectrum of the force field. The Fourier components F k ^ of F(x,t) are defined through 

F(x, t) = J2J2 P ^ e<(fcX_Wt) > ( 2 - 27 ) 

k LO 

with k = n5k, oj = m5u, 5k = 2n/L, and 5u = 2n/T. Here L and T are the spatial and 
temporal periods of the force field. Passing to the continuous case, the Fourier components 
F(k, to) of F(x, t) are defined through 

F(x,t) = J J dkduF(k,Lu)e l(kx - ut) . (2.28) 
Thus the Eulerian correlation function may be written in the form 

C(x, x', t, t') = < F(x, t)F(x', t') > = 




dkdk'dudu' < F{k,u)F{k',u') > e Kkx+k'x'^t-^'t') _ ( 2 29 ) 



For the homogeneous and stationary case considered herein, the right-hand side depends only 
on the differences x — x' and t — t r ; this is possible only if kl = —k and ou' = —u. This in turn 
means that the average < F(k, u) F(k', u') > should be proportional to 5(k + k') 5(u + ou'): 



< F(k, to) F(k', u') > = S(k, to) 5(k + k') 6(u + u') , 



(2.30) 
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which may be taken as a definition of the force spectrum S(k,u). 

The Eulerian correlation function C(x,t) is thus the Fourier transform of the force 
spectrum: 



C(x,t) = J j dkdL)S(k,u)e i ( kx - ut) . (2.31) 
Taking x = and t = 0, we obtain the normalization 

C(0, 0) = < F 2 > = J J dkdcuS(k,u) . (2.32) 

As an example, for the Gaussian correlation function (I2.13p . the corresponding form of 
S(k,u) is 

S(k,u) = (^) <F 2 > e -A^/4. e -^A . (2.33) 

Under the Corrsin approximation, the momentum diffusion coefficient D(p) (Equation [2.23] ) 
can also be expressed in terms of the force spectrum as 

D(p) = J J dkdcoS(k,uj)G(k,uj) , (2.34) 

or, in the discrete CclSC, clS 

D (P) = J2J2 \hA 2 G(k,u) . (2.35) 

k ui 

Here the function (which has the dimension of time) 

G(k,iu)= dt dxP(x,t)e ikx{t) - iuJt = dt < e ikx ®-*" > (2.36) 
Jo J Jo 

is the Fourier representation of the propagator describing the spatial transport of particles, 

i.e., the Green's function of the spatial transport equation; its form depends on the (generally 

stochastic) behavior of x(t). Generalization of the above results to three dimensions is given 

in Appendix A. 



3. Categorizing Acceleration Scenarios 



The models constructed below are all obtained from Equation (I2.23P or equivalently 
from Equation f|2.34j) . They are grouped, in order of complexity of the equation for P(x,t), 
as resonant, non-resonant, or resonant-broadened. 



3.1. Resonant acceleration models 

In resonant acceleration models, the particle trajectory is free-streaming; the stochastic 
equations of motion are p = F(x, t), x = vt. This means that the probability function P(x, t) 
satisfies 

dP(x,t) 8P(x,t) , \ r / \ /o ^ 

-Lll+.-i-li = *(,)*(*), (3.1) 

with solution 

P(x,t) =5(x-vt) . (3.2) 
Equation fl 2 . 2 3 f) then straightforwardly gives 

poo p poo 

D(p)= j dt dxC(x,t)5{x-vt) = / dtC(vt,t) . (3.3) 
Jo J Jo 

In this case the Lagrangian correlation function is related to the Eulerian correlation function 

by Cl(J) = C(vt,t). Note that when v — > 0, Cx(t) — > C(0,t), therefore the Lagrangian 

correlation time tl becomes equal to the Eulerian correlation time r, and D(p) — >■ r < F 2 >. 

In the Fourier representation, Equation ( 12.36!) shows that the propagator is 

POO 

G(k,u) = / dte i{ - kvt ~"V . (3.4) 



o 



As it stands, this integral is not convergent when t — > oo. To evaluate the integral, we 
temporarily add an infinitesimal damping factor v = 0+ in the exponential, so that 



G(k,u)= dte^ kvt -^e- ut = ^ 7r5(cj-kv)+iV — —) , (3.5) 

J [kv — lo) + iv \oj — kv J 

where V is the principal value. The momentum diffusion coefficient (I2.34p takes the form 
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D(p) = 7r / / dk duo S(k,uo) 5{uo — kv) = ir / dkS(k,kv) . (3.6) 



Some insight into the diffusion coefficient D(p) may be earned by normalizing time and 
space coordinates by the Eulerian correlation time and length r and A, respectively. If we 
then define dimensionless frequency and wavenumber by uo = ur and k = Xk, respectively, 
the delta-function term in the propagator (13.51) becomes 

7r 8(co — kv) = r g(k,co] 6) . (3.7) 
Here the dimensionless propagator 

g(k,S5;9) = it5{u-6k) (3.8) 
is a function of the dimensionless momentum variable 

= ^=^, (3.9) 
A/r p 

where the reference momentum p = mX/r has been defined. Note that 9 is also the ratio 
of the Eulerian correlation time r to the transport time scale X/v . The expression (13. 6ft for 
the diffusion coefficient becomes 

£>(p)=(-L) f [ dkduS{k,u)rg(k,u;9) . (3.10) 



Then, noting from Equation (I2.32p that, by definition, 

<F 2 >= (J-^J J J dkduS(k,u) , (3.11) 
we can eliminate the the factor Ar between Equations (I3.10p and (13. lip , yielding 

D{p) = ri{9)<F 2 > , (3.12) 

where the dimensionless function 

£(6) = n d kduS{k,u))g (k,u;6) 
f f dkdu S(k, uj) 
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Since, by Equation (I2.17p . D{p) = tl < F 2 >, Equation (13.121) shows that the Lagrangian 
correlation time is given by tl = r^(6). Hence the function £(#) is simply the ratio of the 
Lagrangian correlation time tl to the Eulerian correlation time r. Further, according to 
Equation (I3.13p . this ratio depends only on the dimensionless momentum variable 9 = p/po 
and on the form of the (dimensionless) force spectrum S(k,oj). 

We illustrate with the Gaussian spectrum (12.331) . which, in terms of dimensionless vari- 
ables, is 

S(k,u) = S e- l2/47T e-* 2/7T . (3.14) 
Substituting this spectrum in Equation (I3.13p . we obtain 

m = JTTleW ' (3 ' 15) 

so that 



D{p) = [i + igSp ' (3 ' 16) 

In the initial stage of the acceleration process, the particle is moving sufficiently slowly to 
feel only the temporal variation of the force: Tl ~ t when p <C po. If < F 2 > is also 
independent of p, D{p) is a constant and so, by Equation (I2.19p . < p 2 > oc t, as in the 
Obukhov model. However, as the particle gains speed, it eventually becomes fast enough 
that it becomes less sensitive to the temporal variation of the force and more to its spatial 
variation. The Lagrangian correlation time becomes of the order of the transport time 
scale which is inversely proportional to the particle velocity, i.e., tl oc p~ l . If < F 2 > is 
independent of p, then we can use Equation (I2.18P and the defining relation r acc = p 2 / D{p) 
for the acceleration time r acc to obtain, in the nonrelativistic regime where the particle energy 
E = p 2 /2m, 



D{p) cc p- 1 ; r acc cx p 3 cx E 3/2 . (3.17) 

It should be noticed that when the force field is the gradient of a potential so that it is the 
potential rather than the force has a Gaussian correlation function (see Appendix A), then 



D{p) cc p 3 ; r acc cx p 5 cx E 5/2 



(3.18) 
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3.1.1. Resonant Acceleration by Waves 

For a wave-acceleration process, we include the dispersion relation u(k) in the Fourier 
spectrum of the force spectrum S(k,u), i.e., 



S(k,uj) = S{k)6[L0-L0{k)} . 
From Equation (13.61) . this leads to a diffusion coefficient of the form 



(3.19) 



D{p) =2tt J dk S{k) 5[u(k) - kv) , 
where the factor of two takes into account waves propagating in both directions. 



(3.20) 



The nature of resonant acceleration by waves is as follows: a particle of momentum p 
interacts with a wave if the particles's velocity is equal to the phase velocity of the wave: 
v = u(k)/k = V p . However, stochastic resonant acceleration can occur even in the absence of 
a dispersion relation. Physically, this occurs because a particle still interacts with a Fourier 
component of the force field such that u/k = v, even though u and k may not be related 
by a dispersion relation uj{k). Thus, particles are stochastically accelerated provided the 
modes which are experienced by the particles along their trajectory have random phases. 
The existence of a broadband spectrum of field fluctuations, in the frame co-moving with 
the particles, with Fourier components not related by a dispersion relation may be identified 
with strong turbulence; this situation relaxes some of the restrictions imposed by a linear 
dispersion relation. 



Sturrockl (119661 ) has shown that (a) by taking into account the dispersion relation 



u = u(k) in the force field spectrum and (b) by enforcing conservation of the total combined 
energy of the field and the accelerated particles, the self-consistent quasil i near diffusion equa- 
tions for wave-particle interaction are recovered (jSagdeev &: Galeev 1969 : Escande fc Elskens 



20031 ). The en ergy gained by t he resonant particles is simply lost by the waves through Lan- 
dau damping. ISturrock! (119661 ) uses as an example an electric force produced by Langmuir 
waves, but the idea readily generalizes to other forms of wave. A detailed discussion of the 
self-consistent quasilinear equations, with e mphasize in the role of resonant and non-resonant 
particles, can be found in iKaufmanl (119721 ). 



In resonant acceleration by plasma waves, the momentum diffusion generally depends 
on the spectral properties of the turbulence. This may be seen by taking a wave frequency 
independent of k, i.e., u = uq, and a power-law spectrum S(k) oc k~ q in (13.201) . then 
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D{p) oc p 1 q which involves the spectral index q of the turbulence^- By Equation ( I2.19p . in 
such a scenario, the particle energy grows as 



Again, it is emphasized that Equation ( 13. 6 p is more general than Equation ( I3.20p . The 
former is valid for an arbitrary correlation function rather than one which describes only 
oscillations and which therefore, for a given k, has a peak at the frequency u(k) of the 
oscillations. 



3.1.2. Remarks 

The domain of validity of the quasilinear approximation is obtained by considering 
the dimensionless Kubo number K = t\J < F 2 >/po that enters the normalized momentum 
equation dp/dt = KF, where the length, time and force have been normalized to r, A and 
V < F 2 >, respectively. When K -C 1, we are in the domain of validity of the quasilinear 
approximation. This inequality limits the typical amplitude of the force field to 



A 



(3-21) 



and gives quantitative meaning to the "weak field limit" discussed by lSturrock! (119661 ) . Also, 
when the spectrum of modes is discrete, their amplitude must be sufficiently large to permit 
resonance overlap, and hence stochasticity, otherwise dynamical trapping of particles become 
substantial and acceleration is inhibited. The problem with spectral discreteness can be 
understood by recalling that the momentum diffusion coefficient in such a case takes the 
form 



D{p) 



dt 



E 



F k ,u,\ 2 exp[i(kv - u)t] 



(3.22) 



which shows that the wave-numb er sum inside th e time integral does not decay with time, 
but instead exhibits recurrences (jKrommesI |2002| ) . These recurrences are a consequence of 
the periodicity of the force felt by the particles along their trajectory in the quasilinear 
approximation x = vt. This effect is similar to (but not equivalent to) the dynamical 
trapping of a particle oscillating in a periodic electric field where both x(t) and p(t) are 
periodic functions of time. The momentum diffusion coefficient can also be written in the 
form 



'This is widely studied in the context of high-frequen cy turbulence of Langmuir waves, which frequency 



const and the spectrum S(k) oc k q (see, e.g.. lTsvtovichlll995l) . 
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D(p) = nJ2 \F k ^\ 2 S[u - kv] . (3.23) 

k,u> 

This form is somewhat is pathological since diffusion occurs only over a set of separated 
points in velocity space. As a result, the distribution function is not affected by the diffusion 
process and acceleration does not occur at all. 

One way in which the singular nature of the diffusion process can be removed consists 
of broadening the resonant propagator tt5(u — kv) into a suitable, say Lorentzian, form: 



OM-J2\^\ \ u _; v)2+1 , 2 ■ (3-24) 



k,U> 

which may be obtained by replacing uj — > u + iv with v > in Equation ( 13.22D . see also 
Equation ( 13. 5ft . Another way of removing the singularity is to simply replace the discrete 
sum in Equation f l3.22|) by a continuous one, i.e., 

Ef dkdu , 

fe,0) 

Identifying S(k,u) = \Fk,J\ 2 /{8k Su), one recovers Equation ( 13. 6ft . 

In the non-resonant or resonance-broadened acceleration models to be presented below, 
the propagator is no longer singular, a consequence of the seed stochasticity that is included 
in the description of the spatial transport of the particles and which is not, in general, 
dynamically related to the accelerating wave field. For this reason, stochastic acceleration 
occurs also when the spectrum of the wave field is singular (discrete or monochromatic) in 
these models. 



3.2. Non-resonant acceleration models 



Non-resonant acceleration models are base d on the approximat i on tha t the spatial trans- 
port is also stochastic; see, e.g., the review by iBykov fc Toptyginl ( 11993b! ) . In the standard 
case, this transport is diffusive, i.e., the particle's position varies according to a Gaussian 
process and the mean square displacement grows linearly with time. In other words, the 
stochastic equations of motion are approximated by p = F(x,t), x = ((t), where ((t) is a 
Gaussian white noise with < C(t)((t') >= 2n5(t — £'), k being a spatial diffusion coefficient. 
The propagator P(x,t) satisfies the diffusion equation 
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°jm- K »™ =S(x)m , ,3,6) 

with the standard Gaussian solution 

P(x,t) = -^L= e~ x2/AKt . (3.27) 

Particles perform a standard random walk, with < x 2 >= 2nt; the characteristic diffusion 
time that is associated with the correlation length A is t d = \ 2 /k. In Fourier space, the 
propagator 



G(k,u) = r dt < e *[M*)-«t] >= / dte -iu* e -#<*>/2 = __L_ = V -±—^- (3.28) 

where the damping factor v& = nk 2 , a parameter associated with the spatial diffusive trans- 
port. For a force field with a continuous spectrum of Fourier components, the diffusion 
coefficient becomes 

D(p) = J J dkduS(k,u) u2 ^ Kk2? , (3-29) 

and involves a Lorentzian of width v^ 1 = (/cfc 2 ) -1 . A similar expression exists for force fields 
that consist of discrete sum of Fourier components: 

If we again normalize time and space by the Eulerian correlation scales r and A, then 

kJc 2 

r 9 (k,u;C) , (3.31) 



co 2 + (nk 2 ) 2 
where the dimensionless propagator 



= (3 - 32) 



is a function of the dimensionless parameter 



(3.33) 
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the ratio of the Eulerian and transport time scales. Again, we can write the momentum 
diffusion coefficient in the form 



D(C) = r£(C)<F 2 > , (3.34) 
where the non-dimensional parameter 

f f dk dcu S(k,u) q(u, k; () 

f (0 = ^ ~ J J ^J ,s; . (3.35) 

f f dkdu S(k, oj) 

When ( < 1, i.e., r <C td, then the propagator reduces to the delta function, g(k,u;( *C 
f) — > tt5(u), as in the Obukhov model. On the other hand, when £ >> 1, i.e., r ^> td, 
the dimensionless propagator becomes g(k,u;( ^> 1) — > ((k 2 )' 1 . Exact forms for may 
be obtained for a prescribed spectrum S(k, a;) of the force field. However, it is sufficient to 
note that when ( <C 1, oc C°, and that when ( > 1, ((() oc Thus, the results are 
summarized as follows: in the weak diffusion limit (r <C r D ), 

D{p) ~ r < F 2 > , (3.36) 
while, in the strong diffusion limit (r ^> td), 

D(p) ~ r D < F 2 > . (3.37) 

As for resonant acceleration, the Lagrangian correlation time interpolates between the Eule- 
rian and the transport time scale, the latter being now the diffusive time scale. Therefore, if 
the force field and the spatial diffusion coefficient k (and hence the transport time scale td) 
are independent of momentum, in both weak and strong diffusion limits D{p) is independent 
of p. 

To summarize, for non-resonant acceleration, in the nonrelativistic limit, 

Dip)^p°; r acc ocp 2 oc £ . (3.38) 

3.2.1. Fractional- Order Fermi acceleration 



Anomalous momentum gain with < p 2 > growing faster than linearly in time can 
occur in non-resonant models, even when the force field and the transport time scale are 
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independent of momentum. This may happen for a a static (d/dt = 0) force field, i.e., 
S(k,u) = S(k) S(uj), for which by Equation (|3.29p . 

j . If S(k)dk 
D(p) = -J . (3.39) 

Moreover, let us take the force field to be scale-free, for instance 



S(k) = So = constant , (3.40) 

so that the integral for D(p) diverges at small k like k~ x and D(p) — >■ oo. It follows that for 
a scale-free force, the momentum gain process is super-diffusive: 



< p 2 > oc ft 3 ; $ > 1 (3.41) 

at all times. However, if the force has a characteristic scale A, the super-diffusive process is 
limited to times less than the time scale = A 2 / k taken for spatial diffusion over a length 
A to occur. Matching the super-diffusive momentum gain regime < p 2 > ~ t" to its diffusive 
counterpart < p 2 > ~ At ~ y/KT^t at t = td, we see that ~ t^ 2 td, so that (3 = 3/2. 
Hence 



< p 2 >~ t 3/2 ; t<r D ; <p 2 >~ t ; t > r D . (3.42) 

Summarizing, 

D{p) -> oo ; r acc oc p 4/3 oc E 2 / 3 . (3.43) 

Superdiffusive momentum transport occurs in scale-free force fields and is the preasymp- 
totic regime of non-resonant acceleration models involving frozen fields. The reason is that 
the random motion of the particles results in multiple returns to any given position. Because 
of the static nature of the force field, the particle experiences the same value of the acceler- 
ating field at each visit to the same position. Such recurrences introduce correlations in the 
stochastic acceleration process that are absent in a non-static environment, rendering the 
situation more akin to a systematic one, and hence faster than diffusive. We may thus call 
such a stochastic acceleration process Fermi acceleration of fractional order. The essential 
difference between second-order Fermi acceleration and fractional-order Fermi acceleration 
is as follows. For the second-order process D(p) exists. It may produce anomalous (superdif- 
fusive) behavior because of the dependence of D{p) on p [e.g., D(p) oc p; < p 2 >oc t 2 ; see 
Equations (I2.18P and ( I2.19j) ] : however, the evolution of f(p,t) is still given by a standard 
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Fokker-Planck equation. However, in fractional-order Fermi acceleration, the anomalous 
superdiffusive behavior in momentum space is due to the fact that the momentum diffu- 
sion coefficient is infinite, and the the equation governing the evolution of f(p,t) is not 
the Fokker-Planck equation. Stochastic acceleration model s base d on the fractional Fokker- 
Planck equation have been discussed by lBian fc Browning! (120081 ). 



3.2.2. Nonresonant Stochastic Acceleration by Turbulent Waves 

When to and k obey a dispersion relation, it follows from writing S(k,u) = S(k) 5[u — 
u)(k)] that the momentum diffusion coefficient, when it exists, is 

/nk 2 
iis «^TW' (3 ' 44) 

Taking the simple dispersion relation u(k) = V p k, and defining the dimensionless parameter 

Cw = ^=^, (3.45) 
with Tw = ^/V p , yields the following form for the momentum diffusion coefficient: 

D(Cw)=t d £(Cw) <F 2 > . (3.46) 
Here the dimensionless function 

= fdkS(kJ g (k-,Cw) 
JdkS(k) 

and 



g(k; Cw) = ,~ ~ • (3.48) 
Cw k2 + kA 

Again, an exact expression for £(Cvk) ma y be obtained for a specified form of the wavenumber 
spectrum S(k). In the limit (w <^ 1, g{kXw < 1) k~ 2 , leading to £,((w "C 1) — > 
constant. In the limit (w 3> 1, g(k,(w 3> 1) — > ( w 2 an< ^ thus £(Cw 3> 1) — » ( w 2 • These 
results are summarized as follows. When t d <C t w , then 



D{p) ~ r D <F 2 > , 



(3.49) 
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while, in limit td 3> Tw, 



In both cases, therefore, 



D {p) ~ — < F 2 > 



(3.50) 



D(p) oc p°; r acc <x.p 2 <x. E 



(3-51) 



3.2.3. Nonresonant stochastic acceleration by a monochromatic wave field 



We have already given - Equation (I3.30P - a general expression for the momentum 
diffusion coefficient which results from non-resonant acceleration by any periodic force field 
consisting of a discrete spectrum of Fourier components. Here we discuss the illustrative 
example of a monochromatic mode with frequency Uq and wave- number fc . From Equa- 
tion (I3.30p . the momentum diffusion coefficient is given in this case by 



D(p) 



I F I 2 
\ r ko,uio I 



w + 



M) 2 



(3.52) 



and the frequency can be taken to obey the linear dispersion relation ujq = V p ko. This result 
provides a simple explanation of the reason why particles are stochastically accelerated by 
adiabatic compressive forces F = — |pV • V (see Section POl) . even if the flow V consists 
of smooth and regular com pressions and e xpans ions. This mechanism was called diffusive 
compression acceleration by iJokipii et all (120031 ). where it was noted that the acceleration 
mechanism has similarities with both second-order Fermi and diffusive shock accelerations; 



see also IZhangl (120101 ). Therefore, a spatially periodic compressive flow produces the same 
kind of stochastic acceleration as a strongly turbulent flow or a shock flow. A legitimate 
question is thus: what is the form of the Lagrangian correlation function, or equivalently, 
what is the shape of the frequency spectrum of the monochromatic wave field in the frame 
of the random walking particle? A simple calculation shows that this spectrum is broad and 
has the Lorentzian form: 



fcoM) I 



71 



(co - u ) 2 + («fcg) 2 ' (3 ' 53) 
Note that this spectrum is similar to the one of an harmonic oscillat or in which t he frequency 
is randomly perturbed in time, i.e., the Anderson-Kubo oscillator ( IKubolll963l ): 
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tp(t) = iu(t)ip ; u(t) = uj + Au(t) . (3.54) 

The frequency perturbations are assumed to have a zero average and a correlation function 
given by 



<Au(t)Au(t')>=—exp(-\t-lf\/T) . 

T 

In the white noise limit r — > 0, the spectrum is indeed a Lorentzian: 



(3.55) 



I{u) 



1 



7T ( W - W o) 2 + ^ ' 

however, in the limit of large r, i.e., (D w t) 1 I 2 ^> 1, we obtain the Gaussian spectrum 



(3.56) 



I{u) = I — — I exp 



-(a; - up) 



(3.57) 



showing how the broadening of the line shape depends on the correlation time of the fre- 
quency perturbations. 



3.2.4. Remarks 

In its original version, the quasilinear diffusion equation is a Fokker-Planck equation that 
describes the diffusion of particles in momentum space due to the action of turbulent waves. 
The quasilinear formalism can be seen in a general perspective by casting it in Hamiltonian 
form. In absence of the the accelerating field, the particle motion is described by an un- 
perturbed Hamiltonian H which is independent of time. This means that the unperturbed 
motion occurs at constant energy and is rectilinear in an infinite homogeneous plasma. The 
total single-particle Hamiltonian is then H + Hi(t) and the change in the particle's energy is 
a consequence of the time- dependent perturbation H\(t) with Hi/H ~ e <^ 1. In standard 
quasilinear theory, the stochasticity comes from Hi(t), i.e., from the assumption of random 
phases in the broad spectrum of the accelerating field. The smallness parameter e expresses 
the condition that the amplitude of the force is small, so that the acceleration efficiency can 
be derived using linear 0(e) theory, hence the terminology "quasilinear." It is not necessary 
to assume an infinite, homogeneous plasma, the fundamental requirement is simply that H 
is integrable. Quasilinear diffusion occurs in action space; the unperturbed motion in angle 
space is also rectilinear, just as that for the motion in coordinate space in a homogeneous 
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medium. The accelerating waves are in resonance with the periodic bounce motion of the 
particles in this case 

However, there is also another form of quasilinear theory which can be applied in the 
complete reverse case to that above, namely when the unperturbed Hamiltonian is non- 
integrable, giving rise to chaotic motion of the particles in the absence of the accelerating 
field. The unperturbed dynamics described by Hq provides the source of randomization for 
stochastic acceleration to occur under the action of the time- dependent perturbation H\. 
Notably, since the unperturbed motion is a seed stochasticity, stochastic acceleration can 
occur also when the perturbation is coherent. Here, we have labeled such quasilinear models 
as "resonant" or "non- resonant." The essence of non-resonant stochastic acceleration is 
therefore simply understood by considering the unperturbed motion to be a spatial random 
walk, with a mean free path of the particles smaller than the correlation length of the field. 



3.3. Resonance-broadened models 



In resonance-broadened models, the spatial transport is a combination of free-streaming 
and spatial dispersion: x — vt + Ax, where Ax is a random walk. Th e idea o f reso nance 



bro adening in the context of acceleration of particl e s was introduced bv iDupred (119661 ) [see 
also lWeinstockl (119681 ) ; lRudakov &: Tsytovichl (ll97ll ); lThomson &: Benfordl (119721 )] . The prin- 
cipal effect of the dispersion can be seen by applying an orbit diffusion < Ax 2 >= 2nt around 
the free-streaming solution; hence, the propagator P(x, t) evolves according to 



dP(x,v,t) dP(x,v,t) d 2 P(x,t) 
dt dx dx 2 

The solution for the propagator P(x,t) is 



5{x)S{t) . 



(3.58) 



P{x,t) 



1 



exp 



(x — vt)' 



(3.59) 



V4nKt V AKt 

When k = 0, P(x,v,t) reduces to 5(x — vt), which is the usual resonance function of quasi- 
linear theory, considered in Section |3~T1 Hence, the essential effect of the diffusive term is to 
produce a broadening of the resonance function. 

Equivalently, the expression 



D(p) 




dk du I dt S(k, u) e~ iuit < e ikx(t) > 



(3.60) 
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is calculated with 



x(t) = vt + Ax{t) , (3.61) 

where Ax(t) is the deviation from free-streaming motion, which is assumed to be a standard 
random walk with < Ax 2 >= 2nt. The momentum diffusion coefficient now involves a 
velocity-shifted Lorentzian: 

f f nk 2 

D(p) = _/ JdkioSik,*) J — W ^ wf . (3.62) 

Scaling as usual, we may write the momentum diffusion coefficient in the form 

D(9;0 = rmO<F 2 > , (3.63) 
where 6 and ( are as defined in Equations (I3.9P and (|3.33p . respectively: 

Q = JJ d kdu S(k, u) g(u, fc; gj Q ^ ^ 
J f dkdu S(k, uj) 

and 



£5; 0; C) = • (3.65) 

(w-^) 2 + (CF) 2 

This result now involves two non-dimensional parameters 9 and (: the model has both a 
resonant and non-resonant character depending on the relative value of these parameters. 
Obviously, when both are small we recover the Obukhov model. 



3.3.1. Dupree theory 



In formal resonant-broadening theories the resonance is broadened by a function R(u — 
kv,i/), not neces sarily a Lorent zian, which tends to 8{u — kv) when the frequency v — > 0. 
In the theory of iDupred (119661 ). the accelerating field is idealized by position- independent 
Gaussian white noise, i.e., the Obukov model, leading to a spatial dispersion with < Ax 2 >= 
(2/3)Dt 3 /m 2 around the f ree-streaming motion. Therefore, the momentum diffusion coeffi- 
cient originally derived by IDupred (119661 ) is 
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D(p) = J dkdco J™ dt S(k, u) e < kvt -^ e - fe2D ' 3 / 3m2 . (3.66) 

Here, the spatial transport in the accelerating field is dynamically correlated with the field 
itself. It is important to notice that a characteristic result of this kind of analysis is an 
implicit, rather than explicit, relation for D, which appears on both sides of Equation (13.661) . 
Formally, this means that in order to obtain the momentum diffusion coefficient D, we need 
to solve an equation of the form 



D(p) = t l {D) < F 2 > , (3.67) 

with tl a function of D; thus D does not straightforwardly scale as the square of the ampli- 
tude of the force. 



3.3.2. Remarks 



To conclude this section, we mention another form of broadening. In the presence 
of wave turbulence, for a given wave-number the frequency spectrum does not have to be 
sharply peaked at the frequency corresponding to the dispersion relation u = u(k). When 
the spectrum of the force field has a peak around co(k) with width rr 1 , the wave spec trum 
can be modeled as a broadened delta-function, for instance (IBykov fc Toptyginlll993bl ) as a 
Lorentzian 



F(k, uj) | 2 = S{k, u) 6[u) - u{k)} -> tt^ 1 S{k) 



-1\2 



(3.68) 



u(kW + (t- 

Broadening of the propagator 5{uj— kv) and broadening of the frequency spectrum 5[u— oj{k)\ 
are formally similar. However, since they involve two different processes they may still be 
assigned two distinctly different physical meanings. The first is a result of the spatial trans- 
port of particles; the second is a spectral property of the turbulent field involving for instance 
non-linear interactions between the waves themselves. In the context of particle transport, 
various mechanisms of reson ance broadening, including w ave d amping or dynamic al turbu- 
lence have been discussed by ISchlickeiser fc Achatzl (119931 ) and iBieber et al.l ( 119941 ) . 
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3.4. Non-Gaussian spatial transport 

Up to now we have assumed that the stochastic part of the spatial transport was 
described by a standard diffusion equation. However spatial transport can described by 
a variety of mode ls, including the Chapman-Kolmogorov and the telegrapher equations 



( IWebb et al.l 120061 ). Let us take, for instance, the propagator to be the solution of the 



fractional diffusion equation 

£>? P{x, t) = K Df xl P(x, t) . (3.69) 
This produces anomalous dispersion^ with 

< x 2 > oc t 2 ^ r (3.70) 

and a corresponding non-Gaussian probability distribution functions P(x,t). For instance, 

resonance broadening by super- diffusive transport (Levy flights) can be modeled by the 
following advection-diffusion equation: 

dP(x.t) dP(x,t) ^ r „, , , n . 

dt~ V dx~ = K ^' • (3 - 71) 
Since the Fourier transform is J-(DTj{x)) =\ k | r f(k), we obtain 

D (P) = / J»toS { k, U ) Hk) _ff +mr)7 - (3.72) 

The corresponding transport time scale is correspondingly different from the diffusive timescale 
v~^ x found earlier. 



3.5. Anisotropic spatial transport 

We now discuss a simple variation of the cases discussed above. We still assume the 
existence of a strong guide field B directed along the x-axis and the existence of a stochastic 
force parallel to the magnetic field F(x, y, z, t), say a parallel electric force or magnetic mirror 
force. The particle transport in this field is taken to be mainly free-streaming along the 
magnetic field but diffusive across the field. Hence, the propagator P(x, y, z, t) satisfies the 
advection-diffusion equation 



6 The standard Gaussian case is recovered with (3 = 1 and T = 2. 
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^^ + ,|,^^ = «xViP(W), (3.73) 

where Vj_ is the two-dimensional Laplacian in the directions perpendicular to the streaming. 
The solution is similar to Equation H3 .591) . with n± being the cross- field spatial diffusion 
coefficient. The corresponding momentum diffusion coefficient is given by 



Din) = 2. Jk ±dk± J i h jtoS(k M {u _^*f- (K ^ , (3.74) 

where the turbulence is also assumed to be isotropic in the plane perpendicular to B . As 
expected, perpendicular diffusion alone can act to broaden the resonance between the modes 
and the particles. Moreover, if these modes are waves with frequency u which is related to 
the wavenumber k through the anisotropic dispersion relation u = u(kj_, kn), then the force 
field spectrum S(k±, k\\, co) = S(k±,k\\)S[cu — cj(k±,k\\)], and 



K\_k\ 

[u(k±, k\\) - k\\v\\} 2 + (n±k] 



D(p\\) = 2tt / k± dk L I dk\\ S(k ± , fey) LJU U A ,_^ 12 ; ^ u2 . (3.75) 



4. Forms of the Accelerating Force Field 

We now move on to consider various possible type of force fields and the resulting forms 
of D(p). Recall from Equation (I2.17P that in general D(p) = tl < F 2 >, and it should 
be noted that the dependence of D on p can arise both through the Lagrangian correlation 
timescale tl and/or because < F 2 > is itself momentum- dependent. 



4.1. Drift-kinetic equations — magnetic mirror and parallel electric forces 

Since turbulent energy is generally a decreasing function of frequency in magnetized 
plasmas, it is reasonable to assume that low-frequency turbulence is most efficient at accel- 
erating particles. For electromagnetic fluctuations that satisfy uj Ud, where u C i is the ion 
gyrofrequency, the motion of particle gyrocenters obeys the drift kinetic equations 

dX , fi . 

— = D||b + v ExB , (4.1) 
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dv\\ q 

-JL = ^-^V„B, (4.2) 

= . (4.3) 

Here v ExB = c(E x B)/S 2 is the E x B drift velocity, b = ~B/B, B —\ B | and the magnetic 
moment 

= f| • (4-4) 

The corresponding kinetic equation is thus 

^ + hb + v ExB ) -V/ + (^|| - /, W V||S) = . (4.5) 

This shows that low-frequency turbulence with a; <C w C i accelerates particles both through 
the fluctuating parallel electric force 



F|| = qE\\ (4.6) 
and the fluctuating magnetic-mirror force 

*ii = ljF v '' B ■ (4 - 7) 

These two parallel forces are responsible for Landau damping (or transit-time damping) of 
various plasma waves via resonant acceleration of the particles. In general, however, particle 
acceleration through these stochastic forces can be either resonant or non-resonant. In the 
latter case, zero- frequency modes are absorbed through particle acceleration by a process akin 
to Ohmic or viscous damping. Acceleration by these forces is also intrinsically anisotropic. To 
produce an isotropic distribution function, therefore, a high level of pitch-angle scattering 
may also be required. Thus, a scenario can be envisaged where one component of the 
turbulent fluctuations produces, on a fast time-scale, the isotropization of the distribution 
function via pitch-angle scattering. This results in a diffusive transport of the particles at 
constant energy with a spatial transport coefficient which is in general a function of the 
particle speed. Then, another component of the turbulence acts on this diffusive motion 
to produce the energy chang e, resulting in stochastic acceleration of the particles. This is 



essentially the iParkerl ( 119651 ) formalism, where energy changes are a result of large scale 
compressions. 
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4.2. Finite Larmor-radius effects 

In the drift-kinetic approximation, the electromagnetic field fluctuations are an average 
of the true fluctuations over the fast time scales associated with the gyromotion of the 
particles. Therefore, in this approximation, the particles feel an effective low-frequency field 
which acts on their gyrocenter. For k±p <C 1 (where p = v±/ou c is the gyroradius), the 
particle position and the gyrocenter position are essentially the same. However, for k±p > 1, 
the effect of the finite gyroradius must be taken into account. For circular motion, the 
gyrocenter position X is related to the particle position x via 

X = x-p( ¥ >), (4-8) 
where f is the phase and p is the gyroradius vector defined by 

b ° x v ^ (a n\ 

P = • (4-9) 

Upon introducing the Fourier representation F(x, t) = J J dkdu F(k,u) e J ( k - x -<^*) of an ar- 
bitrary field F(x, t), one finds that 



<F(x,t) > (p =F(X,t) = J (ik(fwF(k 1 u)e i(k ' x - ul) ^- J dip . (4.10) 

Since k • p((p) = k±p cos ip, one can make use of the generating function for Bessel functions 

e iz cos v = Y^t-ovi™ J n (z) e inv to obtain 

<F(x,t)> v = J dkdu J (k ± p)F(k,u)e l(l1 ^-^ . (4.11) 

We see that the effective field felt by the particle is reduced by a factor Jo(kj_p) which 
depends on the perpendicular wave-number; this correction is a finite-Larmor-radius (FLR) 
effect. The FLR effect is most simply accounted for in the Fourier representation where it 
enters as a simple renormalization of the force field spectrum 

S(k, u) -)• 4 {k ± p) S{k, oo) . (4.12) 
Moreover, we can define the gyrocenter propagator as 



/•CO 

G(k, u) = / dt < e lk X «-^ > 
Jo 



(4.13) 
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The FLR modified momentum diffusion coefficient takes the form 



D = J J dkdu J$(k ±P ) S(k,u) G(k,w) . 



(4.14) 



In the case where the unperturbed guiding-center motion is rectilinear at constant speed vn 
along the magnetic field, k • X(t) = k\\V\\t, and therefore 



D (PbP±) 



7T 




dk du Jn 



S(k, uj) 5{to — k\\v II 



(4.15) 



The momentum diffusion coefficient (I4.15P involved in the Landau resonance between the 
field and the particles is written here in terms of the parallel and perpendicular momentum 
components p\\ and p±. Alternatively, it can be expressed in terms of the two independent 

'pj_ + pfi is the magnitude of the momentum and fx = p\\/p 



variables p and fx, where p - 
is the direction pitch-angle cosine. For computational applications, it can also be us eful to 
replace the Bessel function by its Pade approximation as (see, e.g., iPress et al.lll992l ): 



J (k ± p) 



1 + (k ±P ) 2 /4: 



(4.16) 



It is important to include the FLR renormalization of the field spectrum when describing 
the acceleration and transport of particles by low-frequency turbulence in situations where 
the field fluctuations have significant en ergy at perpend i cular wave numbers, i.e., at wave- 
lengths of the order of or larger than p. iBian fc Kontarl (120101 ) have discussed the spectral 
structure of the parallel ele ctric field associated with strong anisotropic Alfvenic turbulence 
(IGoldreich fc Sridharlll995l ). showing that the latter attains maximum energy when kj_pi ~ 1. 
The frequency of s uch anisotropic flu c tuations remains well belo w the ion cyclotron frequency 
even at k±pi ~ 1 (IHowes et al.ll2006l ; ISchekochihin et al.l 120091 ). At the cyclotron frequency, 
another FLR effect comes into play which is the interaction between the gyromotion of the 
particles and perpendicular electric fields. The full expression for the quasilinear momentum 
diffusion tensor, including bo t h the Landau resonan c e and the gyroresonances can be found i n 
Kennel fc Engelmannl (119661 ) , iJaekel &: Schlickeiserl (119921 ) and ISchlickeiser fc Achatzl ( 119931 ) . 
New forms of the gyropha se averaged Fokker -Plan ck coefficients for particl e trans port have 
recently been discussed by ISchlickeiser! (120111 ) and lCasanova fc Schlickeiserl (]2012f ) . 



-32 - 



4.3. Acceleration by magnetic mirrors produced by fast-mode waves 

A popular model for flares is the acceleration by the fluctuating magnetic mirror force 



assoc i ated with fast mode turbulence ( Kulsrud fc Ferr arill 1 9 7 It lAchterberg 



1996, 1997; Schlickeiser fc Miller 1998 



1981: Miller et al. 



Yan et al. l2008h . For resonant acceleration by waves 



which possess a parallel electric or magnetic mirror force (i.e., waves that are subject to 
Landau or transit-time damping), we have 



4vr z 



k±dk± / dk\\S(k±,ku) 5[u(k±,k\\) — k\\v\\] 



(4.17) 



where we have neglected the FLR effect discussed above. 



For such a resonant acceleration process the force [Equation ( 14. 7p ] is the gradient of 
a potential. Hence, by Equation (I3.18p . r L oc 1/v 3 . Since the turbulence is weak, the 
Lagrangian time scale should also be proportional to the Alfven time ta = A/Va, where A 
is the characteristic length scale, so that Tl ~ (Va/v) 3 (A/Va). Moreover, the magnitude of 
the force (14.71) is F ~ m e v 2 (1/A) (B/B ). Hence we expect that 



D(p)=r L F 2 



m,v;i i ■ : i — 



Br 



P 



(4.18) 



The exact expression may be obtained as follows. Using the expression (14.71) for the 
magnetic mirror force and the dispersion relation u> = kVA for fast-mode waves (where Va is 
the Alfven speed), we obtain 



ir 2 ml 4 



B 2 



j kldk ± J ' dk\\kfS B \\(kj_,k\\)6(kV A -k\\v\\) , (4.19) 

with Sb,, {k±, k\\) the spectrum of the field strength fluctuations B\\. This diffusion coefficient 
may also, of course, be written in terms of the two independent variables p and fi. As- 
suming that efficient pitch-angle scattering maintains near-isotropy of the electron velocity 
distribution function, the momentum diffusion coefficient may then be averaged over [i: 



D(p) 



d(j,D((jt,p) 



(4.20) 



with the result (IMiller et al.lll996l . Il997l ) 



D(p) 



1 ~ v 21 

e VA 



16 



B 



A \B, 



P 



(4.21) 
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consistent with the estimate (I4.18p . Note that for this mechanism, 



D(p) oc p; r acc oc p oc E 1/2 



(4.22) 



Transit-time acceleration by fast modes was als o studied by ISchlickeiser &: Miller! (Il998[ ) 
who use the full expression for the diffusion tensor (IJaekel &: Schlickeiserlll992l ). They show 
that for flat (1 < q < 2) turbulence power spectra, the acceleration time scale increases as 
r acc oc E( 3 ~~ q ^ 2 , while for steep spectra (2 < q < 6) they recover r a f . r oc E 1 / 2 at nonrelativistic 



energies, see Equations (102) and (103) in ISchlickeiser & Millerl (11998f ). lYan et all (120081 ) 
also discuss transit-time acceleration by fast modes in the low-/3 plasma of the solar corona, 
they include a high-A; cutoff scale owing to collisionless damping of the fast waves resulting 
in r acc oc E 1/2 . 



4.4. The Parker equation — adiabatic compression 



The intera ction of parti cles with low-frequency compressible turbulence can also be 
modeled by the iParkerl ( 119651 ) equation, in which particles are accelerated by a compressive 
force. This force or the rate of particle momentum change follows an equation which is 
obtained from an adiabatic thermodynamic process 



V- Vp 



dp 1 
~dl ~ ~3 

The evolution of the isotropic part of the distribution function is given by 



(4.23) 



df 1 df 

^ + v-v/ = v- K v/ + - P v-v^ 



(4.24) 



Notice that spatial diffusion forms the basis of the Parker equation, such that in this frame 
work stochastic acce l eratio n is always non- resonant (IBykov fc Toptygin!ll993bl ; IZhang fc Lee 
20 111 ; iJokipii fc Led 120101 ) . When the Peclet number is small (VL/k 1) and when 
V <^ ^Jk/t, turbulent velocity fluctuations are unimportant for spatial transport, while 
still producing acceleration. However, in general the presence of turbulent velocity fluc- 
tuations V enhances the spatial transport coefficient and, hence, the turbulent transport 
coefficient K t is given by 



n t = k + r < V 2 > 



(4.25) 
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The momentum d i ffusion coefficient in the case of weak turbulence was written in, e.g., 
Bvkov fc Toptveinl jl993bh as 



(4.26) 



with J dkS(k) =< (V.V) 2 > . [This result can be obtained by blending Equations (I3.29P 
and ( jXSSD -] 

When the correlation time r r — > oo, the momentum diffusion coefficient becomes 



D(p) 



V 
9~ 



dkS(k) 



K t k 2 



(4.27) 



[oo 2 (k) + (K t k 2 ) 2 } ' 

The spectral density S(k) is in principle related to the spectrum of plasma density fluctua- 
tions, while the dispersion relation w(k) may correspond to any low -frequency compressiv e 



201 



modes, including an anisotropic spectrum of kinetic Alfven wave s (IBian Sz Kontarl 
Both fast and slow magnetosonic wave modes have been considered ( Ptuskin 19881 ; Chandran 
20031 : IChandran fc MaroiJliooS Icho fc Lazarianil2006h . 



4-4- 1- Example 1: Non-resonant acceleration by adiabatic compressions associated with 

fast-mode waves 

In this case the dispersion relation co (k) = kVA- Assuming isotropic turbulence S(k) dk - 
47ik 2 S(k) dk, it is found that in the strong diffusion limit Kt/X ^> Va, 

D{p) . ( K t 



A 

while in the small diffusion limit Kt/\ <C Va 



l>2 ^(^) <(V-V) 2 > , (1.2S) 



(- ) <(V-V) 2 > . (4.29) 



Thus, for this mechanism, we obtain the relations 



D(p) oc p 2 ; r acc oc constant . (4.30) 



Note that iPtuskinl (119881 ) has discussed the case of a power-law spectrum, i.e., S(k) oc k q , 
and has concluded that that the scaling of D with p is insensitive to the value of the spectral 
index q. 
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4-4-2- Example 2: Non-resonant acceleration by adiabatic compressions associated with 

strong turbulence 



For the case of strong turbulence (IBykov fc Fleishman! |2009bJ) , 



D(p) = — duj dkS(k, oj) t-= ; — — r-^r 



In the weak diffusion limit (nt *C A 2 /r) 



^Ur<(V-V) 2 > 



while in the strong diffusion limit {n t ^> A 2 /r), 



(4.3i; 



(4.32) 



v V K t . 

For both cases, therefore, we obtain the relationships 



(4.33) 



D(p) oc p 2 ; r acc oc constant . 



(4.34) 



Most non-resonant acceleration models are based on the adiabatic compression relation- 
ship (14 .23 p . Note, however, that compressive waves can non- resonantly accelerate particles 
through both the parallel electric force and the magnetic mirror force. 



5. Energy Spectra 

In the examples discussed above, the momentum diffusion coefficient is a power law, 
i.e., D(p) = Dop a . Therefore, in three dimensions, the distribution function f(p,t) obeys 
the diffusion equation 



df(p,t) _ 1 d 
dt p 2 dp 



D p 



a+2 



df(p,t) 
dp 



(5.1) 



with the normalization provided by J °° 4-7rp 2 f{p,t) dp = 1. Taking the initial condition as 
f(P: 0) = 8(p), the general solution is the stretched exponential 



f(p,t) 



C_ 

47T 



t 3/(2-a) 



exp 



-P 



2-a 



(2 - a) 2 D t 



(5.2) 
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with the constant C = (2 - a )-(i+<*)/(2-°) £> Q 3 /( 2 Q ) r[3/(2 - a)]. This solution is invalid for 
a > 2 since / becomes a monotonically increasing function of p. For arbitrary a, a general 
solution of (15.11) can be found with initial condition f (p, 0) = <5(jo — k), where is the initial 
mom entum. In units such that Pi = 1, we obtain (IBorovsky fc Eilekl Il986l ; iJokipii fc Lee 



2010) 



l<a<2:/(p,t) 



P 



-(l+a)/2 



(2 - a)47rD ^ 
where / = (1 + a)/(2 — a), and 



2p 



(2-a)/2 



(2 -a) 2 A)* 



exp 



1 +p 



2-a 



(2 -a) 2 D t 



(5.3) 



a>2:/(p,t) 



P' 



-(l+a)/2 



2p 



(2-a)/2 



exp 



1 + p 



2-a 



fa -2) 2 D n t 



(5.4) 



(a - 2)47rD t [ — 2 ) 2 ^o^. 
where I — (ot+ l)/(at — 2). In the limit t — >• oo, the a > 2 solution approaches the power-law 
form 



a > 2 : /(p, t oo) oc 



2a-l 



P 



-(1+a) 



(5.5) 



The case a = 2 is special ( IKardashevlll962[ ): 



^•^4^^ /2 ^(-i B0i ) eXP (-S) ■ <5 ' 6) 
These solutions generalize to the case where D is a function of time by replacing the term 
D^t by J* Do(t') dt' in the previous expressions. Thus, power-law spectra can be obtained 
only for a > 2, in which case / ~ the case a = 2 leads to / ~ p~ 3 / 2 . 

For continuous injection of particles in the acceleration region, i.e., when a source term 
5{p — 1) is added to the right side of the diffusion Equation (15. ip . a stationary power-law 
solution may nevertheless result. For such a stationary solution to exist, the momentum flux 
Dop a+2 df j dp must be constant, which requires that 



f(p)~p-(i+«); a >-l . (5.7) 

Note that this has the same power-law spectral dependence on momentum as the (time- 
dependent) a > 2 large-time limit solution (15. 5p . 
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6. Summary 

To summarize, in the stochastic acceleration models discussed above the momentum 
diffusion coefficient D(p) may be written in the equivalent forms 



D(p)= / dtdxC(x,t)P(x,t); D(p) = / dkdcu S(k,u) G(k,u) , (6.1 





where C(x,t) is the Eulerian correlation function of a force F(x,t), S(k,u) is the Fourier 
transform of C(x,t), P(x,t) is the solution of the spatial transport equation, and G(k,u) is 
the Fourier transform of the associated Green's function. Given these general expressions, 
stochastic acceleration models are grouped as resonant, non-resonant or resonance broad- 
ened, depending on the form of P(x,t) or, equivalently G(k,u), which both describe the 
spatial transport of the particles in the acceleration region. For each case, we now summa- 
rize the pertinent equation for P(x,t), its solution, and the resulting G(k,u). 



• Resonant 



dP(x,t) dP(x,t) n 

J +v a =o ; 

at ox 

P(x,t) = S(x-vt) ; D(p) = n dkduS{k,u)5{u-kv) . (6.2) 



• Non-resonant 

dP(x,t) _ d 2 P(x,t) 
dt dx 2 

1 f f Kk 2 

^• t) = ^ "^i/"*" 1 ^' (6 - 3) 

• Resonant-broadened 

This combines elements of both the resonant and non-resonant cases: 

dP(x,v,t) dP(x,v,t) _ d 2 P(x,t) 
dt dx dx 2 



(6.4) 
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Finite Larmor-radius effects can be accounted for by writing the momentum diffusion 
coefficient as 



D(p) = J J dkdu J 2 (k ± p)S(k,u)G(k,u) ; G(k,u) = J dt < e *-X(*)-<arf > ? 

(6.5) 

where G(k, uj) is now the propagator associated with the guiding-center X(i) of the particles. 

The momentum dependence of D(p) arises either because the accelerating force is itself 
momentum dependent (S(k, u>) depends on p) or because the spatial transport is momentum 
dependent (G(k,u>) depends on p). A typical dependence D(p) is a power law 

D(p) = D p a . (6.6) 
Then, since d < p 2 > /dt = 2D(p) = 2D p a , it follows that 

< p 2 >oc t 2 ^ , (6.7) 

for a < 2. For a > 2, < p 2 > can attain infinite values within a finite period of time. Power- 
law distribution functions, solutions of the momentum diffusion equation are obtained only 
for a > 2, in which case / ~ p~( 1+a \ When a < 2 the solutions are stretched exponentials. 
We also noticed that it may occur that the momentum diffusion coefficient is not existing 
because the integral involved in its expression is non-definite. We gave a specific example 
and dumb such acceleration mechanism fractional-order Fermi process. 

All the above models (resonant/broadened/non-resonant) can apply equally well to de- 
scribe stochastic acceleration in a force field which is fragmented, or distributed in clumps, 
as is already known for the resonant interaction between particles and Langmuir waves 



( [Melrose Sz Cramerlll989[ ). Noteworthy, in both resonant-broadened and non- resonant accel- 



eration the spectrum can also be taken to be discrete and even monochromatic. 

In general, the wavenumber k and frequency u are a priori independent quantities. 
However, in the case of acceleration by waves (e.g., MHD waves, plasma waves), these 
quantities are related by a dispersion relation u = u(k). In such a case, 



S(k,u) = S(k) S[u - u(k)] . (6.8) 

For strong turbulence, the frequency spectrum of the force field is broadened around cj(k), 
so that 
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where r c is a turbulence correlation time scale. As remarked above, it is important to note 
that while broadening of the function G(k,u) arises as a transport effect, broadening of the 
function S(k,u) does not. However, the formal consequence is the same - the introduction 
of an additional time scale in the stochastic acceleration model, either as a particle transport 
time scale v^ x = [nk 2 )^ 1 or as a turbulence correlation time scale r c . 

As a rule of thumb, the momentum diffusion coefficient takes the form 



Dip) = r L <F 2 > , (6.10) 

where the Lagrangian correlation time tl is a statistical property of the field in the frame 
comoving with the particles. In all cases, the acceleration efficiency, i.e., < p 2 > as a function 
of time, is estimated using quasilinear theory. Within quasilinear theory, the accelerating 
field acts as a perturbation on the particle motion which occurs at constant energy and is 
either integrable for resonant models or chaotic and diffusive for non-resonant models. 

This general characterization of stochastic acceleration models allows us to use the 
dependence of the diffusion coefficient D on the particle momentum p to compare the ac- 
celeration times, acceleration rates, and energy spectra for a variety of acceleration models. 
Specific expressions for the momentum dependence of the diffusion constant D(p) and for 
the corresponding energy dependence of the acceleration time r acc (E) have been provided 
throughout the text. These predicted scalings can be compared with quantities deduced 
from observation, and so to assess the viability of a given model. Applicable observations 
include the rise time of the radiation field (e.g., hard X-ray bremsstrahlung, gamma-rays) 
produced by the accelerated particles. Hard X-ray measurements al so provide informat ion 



on the specific acceleration rate v sp ecific (particles s -1 per particle; lEmslie et al.ll2008l ). a 
quantity which measures the fraction of the ambient particle population that suffers accel- 
eration to energies equal to or greater than a prescribed energy. With a framework now in 
place to determine the diffusion coefficient for a general stochastic acceleration mechanism, 
the specific acceleration rate can be readily assessed and compared with values deduced from 
analysis of hard X-ray observations. This will be the subject of a future paper. 



A. Generalization of Equations of Motion to Three Dimensions 



Applications of the results of Section [2] are straightforwardly generalizable to the fol- 
lowing equations of motion in three dimensions: 
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p = F(x, t) , mx = p . 
The evolution of the distribution function is given by a Fokker-Planck equation 



(Al) 



df(x,P,t) p , = g 

<9t m ' ' 



Aj(p; 



d/(x,p,t) 
9Pj 



(A2) 



The expression for the momentum diffusion tensor is similar to Equation (12. 23ft and becomes 



Ai(p) 



dt y dxC ii (x,t)P(x,*) 



(A3) 



where the Eulerian correlation tensor of the force field is defined by 



C ij {x,t)=<F i {0,0)F j {x,t) > 
while Equation (I2.34p becomes 



Aj(p) = J J dkdco Sij(k,tu) G(k, 



(A4) 



(A5) 



with definitions of Sij(k,u) and G(k,u) similar to those in Equations (I2.30P and (I2.36p . A 
general decomposition of the Eulerian correlation tensor, assumed to be homogeneous and 
isotropic, is Cij(x,t) = A(x,t)8ij + B(x,t) X{Xj + G(x, t) e^j xi, where 5y is the Kronecker 
delta and e^i is the permutation tensor. The symmetric part of this tensor is invariant relative 
to rotations and reflections. The antisymmetric part changes its sign under reflection and 
therefore is non-zero only for fields possessing helicity, i.e., when F • V x F ^ 0. In a 
similar way, the general decomposition of the spectral tensor of the force field is Sij(k,t) = 
A(k, uj) dij + B(k, oj) kikj/k 2 + G(k, u) e^i ki/k 2 . The particular case of a potential field F = 
— V0, with k x F(k, u) = (longitudinal field), yields 



kk- 

Sij(k,u) = S L (k,u) -j-^- . 
The case of a solenoidal field V • F = 0, with k • F = (transverse field), reads 



(A6) 



Sij(k } uj) = S T (k,uj) 5ij 



(A7) 



Specific geometrical features of the field can be i mportant in studying m omentum and angular 
momentum injections into the particle motion (IChechkin et al.lll994[ ). Here we only focus 
on the energy transfer from the field to the particles. 
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For example, these results can be used to calculate the momentum diffusion tensor 
resulting from resonant acceleration by a the force field which is the gradient of a potential 
F = — V0. In this case, Cij(x,t) = —d 2 C^(x,t)/dxidxj and S^k, cj) = k 2 S<p{k,oj). We 
also assume that (f) has a Gaussian correlation function. Then a calculation similar to that 
leading to Equation (13.161) yields the following expression for the diffusion tensor 



Ai(p) 



l 



PiPj 



Dn 



+ 



1 



( p / p0 )2]3/2 p 2 [l+(p/po)2]l/ 



6y - 



P 2 



(A8) 



with 



r < 



> /A 2 . The isotropic part of the distribution function obeys the diffusion 



equation (in d dimensions): 



df( P ,t) 
dt 



P 



]_d_ 



dp 



(A9) 



where 



Dip) 



D 



;i+p 2 /po) 3/2 

Notice that D(j>) ~ p~ 3 when p > po- This D{p) ~ p 
[Equation (13. 17ft ], is a fundamental consequence of the force being the gradient of a potential 



(A10) 

~ 3 behavior, instead of D(p) ~ p^ 1 
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